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Abstract 

The geometric and electronic properties of Bi-adsorbed monolayer graphene, en¬ 
riched by the strong effect of substrate, are investigated by first-principles calcu¬ 
lations. The six-layered substrate, corrugated buffer layer, and slightly deformed 
monolayer graphene are all simulated. Adatom arrangements are thoroughly stud¬ 
ied by analyzing the ground-state energies, bismuth adsorption energies, and Bi-Bi 
interaction energies of different adatom heights, inter-adatom distance, adsorption 
sites, and hexagonal positions. A hexagonal array of Bi atoms is dominated by the 
interactions between the buffer layer and the monolayer graphene. An increase in 
temperature can overcome a ~ 50 meV energy barrier and induce triangular and 
rectangular nanoclusters. The most stable and metastable structures agree with the 
scanning tunneling microscopy measurements. The density of states exhibits a finite 
value at the Fermi level, a dip at ~ —0.2 eV, and a peak at ~ —0.6 eV, as observed 
in the experimental measurements of the tunneling conductance. 
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1. Introduction 


Two-dimensional graphene has been a main-stream material in both fundamental re¬ 
search and device applications since 2004, [7M mainly owing to its remarkable physical, 
chemical and material properties. [MS] Graphene possesses a unique carrier mobility, mu 
thermal conductivity, mechanical strength, [5112] a quantum Hall effect, [TUlITT] and su¬ 
perfluidity. [l2l[T3] Monolayer graphene is a zero-gap semiconductor since the Dirac points 
have a vanishing density of states (DOS) at the Fermi level {Ep). The electronic proper¬ 
ties can be tuned by changing the doping elements, [?,|29l[30l[33l[3ll[37| geometric curva¬ 
tures. |19p20] laver numbers. |2ip22j stacking conhgnrations. [?J22| mechanical strain, [211125] 
and the applied electric and magnetic helds. [26l[2^ Adatom-adsorbed graphenes can be 
successfully produced in laboratory, such as wide-gap hydrogenated graphenes PI and 
metallic lithium-doped graphenes with a high free carrier density. [?] They are critically 
important in applications of nanoelectronic, nanophotonic, and energy storage devices. 
This work shows how the energetically favorable distribution conhguration of Bi-adatoms 
can be determined by the critical interactions among the substrate, buffer layer, monolayer 
graphene, and adatoms. 

The 2p^ orbitals of carbon atoms are perpendicular to the graphene surface, providing 
a potential chemical environment for doping with various elements with the strong carbon- 
adatom bondings. Experimental research studies have successfully modulated the essential 
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properties of graphene by using dopants, such as O, [221121] H, [301131] N, [321133] transition 


metals, [331135] and heavy metals. [361132] These doping elements induce the novel physical 
properties, including conhguration transformations, semiconductor-metal transition, band- 
gap tuning, adatom-related DOS, creation of magnetic moments, and modulated spintron- 
ics. On the other hand, previous theoretical studies on adatom adsorptions revealed that 
the electronic properties are altered dramatically under various concentrations and distri¬ 
butions. [381 - 130] Moreover, uniform or aggregated adatom distributions play an important 
role in tuning the electronic properties even at the same concentration. Ha Apparently, 
the special distribution of adatoms becomes a critical issue in fabricating high-quality na- 
noelectronic devices. 

To date, bismuth-related systems are some of the most widely studied materials, since 
the dimensionality can enrich the fundamental properties. Bulk bismuth with rhombohe- 
dral symmetry is a semimetal with a long Fermi wavelength and small effective electron 
mass. [32] Its surface states belong to a Dirac fermion gas. SB The three bismuth surfaces: 
Bi (111), Bi (100), and Bi (110) have a higher free carrier density at Ep = 0 compared to 
those of the bulk system. [32] The Bi thin him, in which the surface does not have strong 
chemical reactions with O 2 , appears to be stable up to about 600 K. [33] Furthermore, a 
ID Bismuth nanowire displays narrow band gaps due to signihcant quantum conhnement 
effects. [33] Moreover, bismuth has been explored in the helds of environmental engineering, 
biochemistry, and energy engineering. For example, Bi-based nanoelectrode arrays are used 
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to detect heavy metals, [IS] a polycrystalline bismuth oxide him can serve as biosensor [15] . 


and bismuth oxide on nickel foam covered with thin carbon layers is used as the anode in 
lithium battery. m 

Recently, bismuth adatoms on graphene formed on a 4H-SiC (0001) substrate were 
clearly observed at room temperature. [50l|5l] The corrugated substrate and monolayer 
graphene layers have been identihed by scanning tunneling microscopy (STM). Specihcally, 
a large-scale hexagonal array of Bi atoms was revealed at room temperature. Such adatoms 
are aggregated in triangular and rectangular nanoclusters of uniform size under the further 
annealing process. The scan tunneling spectroscopy (STS) measurements of the dl/dV 
spectrum conhrm the Dirac point, free conduction electrons, and Bi-related structures. 
These results shed light on controlling the nucleation of adatoms and subsequent growth of 
nanostructures on graphene surfaces. On the other hand, there are some theoretical studies 
on the geometric structures and energy bands of Bi-adsorbed and Bi-intercalated graphenes. 
[1HII19] The former have been performed on monolayer graphene without simulation of the 
substrate and the buffer graphene layer, and thus the deformed graphene surface structure 
may be unreliable. [48] The latter is calculated for Bi and/or Sb as a buffer layer above 
the four-layer SiC substrate, which results in an energetically unfavorable enviroment for 
the metal atoms to be adsorbed on the graphene sheet. [49] A comprehensive study on the 
critical roles played by the conhguration of the Bi adatoms, buffer layer, and substrate is 
absent to date. 
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How the energetically favorable distribution configuration of Bi adatoms can be deter¬ 


mined by the substrate effect and the buffer layer is explored in detail by the first-principles 
calculations. This work shows that the adsorption energy strongly depends on the various 
atomic sites and hexagonal positions. The former include hollow, top, and bridge sites, 
and the latter are closely related to the non-uniform Van der Walls interactions. Two 
kinds of different bismuth distributions could be obtained, namely uniform or aggregated 
adatom arrangements with critical energy barriers induced by the buffer layer. The ag¬ 
gregated adatom arrangements are associated with the large-scale hexagonal symmetry, 
Bi coverage, and relatively few vacancies. The optimized geometric structures were val¬ 
idated by the STM measurements. (SHIEI] Hopefully, these rich fundamental features in 
bismuth structures can promote potential applications for use in electronic devices and 
energy materials. Moreover, the main effects of the Bi-adsorption on the density of states 
(DOS) could be understood through a exhaustive comparison with tunneling conductance 
measurements. [sniEi] 

2. Methods 

The first-principles calculations are based on the density functional theory (DFT) im¬ 
plemented by the Vienna ab initio simulation package. [52] The generalized gradient ap¬ 
proximation, within the Perdew-Burke-Ernzerhof functional, [53] is applied to describe the 
exchange-correlation energy of interacting electrons. The projector augmented wave is uti- 
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lized to characterize the electron-ion interactions. The Van der Waals force is included in 


the calculations using the semiemprical DFT-D2 correction of Grimme to correctly rep¬ 
resent the atomic interactions between graphene layers. [51] A vacuum space of 15 A is 
inserted between periodic images to avoid their interactions. The cutoff energies of the 
wave function expanded by plane waves are chosen to be 400 eV. For the calculations of 
the electronic properties and the optimal geometric structures, the hrst Brillouin zones are 
sampled by 3 x 3 x 1 A:-points via the Monkhorst-Pack scheme. The convergence of the 
Helmann-Feymann force is set to be 0.01 eV A 

3. Results and discussion 

First, in order to accurately simulate the bismuth-adsorbed monolayer graphene, the six- 
layer Si-terminated 4H-SiC (0001) substrate is taken into account. The optimized results 
show that the four-layer substrate (region I in Fig. 1), which is reduced from the six-layer 
one, has almost identical geometric properties to the six-layered one. Second, a buffer 
layer, as shown in region II of Fig. 1, is in a periodic ripple shape after relaxation, in 
which the troughs bond with the silicon atoms of the substrate. This periodic corrugation 
is almost identical to that revealed by the STM measurements. [50] Third, the monolayer 
graphene is nearly flat with a slightly extended C-C bond length of 1.50 A, as indicated 
in region III. The interlayer distance between the monolayer graphene and the buffer layer 
varies from 3.21 A to 5.45 A at the crests and troughs, respectively. This difference clearly 
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illustrates the non-uniform Van der Waals interactions between them, and thus dominates 


the distribution of the Bi adatoms. Finally, bismuth atoms can be adsorbed on monolayer 
graphene in self-consistent calculations (region IV in Fig. 1). Up to now, based on the 
different experimental environments, two kinds of adatom distributions have been observed, 
namely a uniform hexagonal distribution and bismuth nanoclusters; these two distributions 
are depended on the adsorption energies. 

The adsorption energy AE, characterizing the reduced energy due to the bismuth- 
adsorption on graphene, is very useful for understanding the optimized geometric structure. 
It is dehned as 

AE = Erotal — Eora — Esuf — Ebi, (1) 

where Erotai, Ecra, Ebu/, and Ebi are the total energies of the composite system, pris¬ 
tine monolayer graphene, buffer layer, and isolated bismuth atoms, respectively. Three 
adsorption sites with higher geometric symmetry are investigated, i.e., the hollow, bridge, 
and top sites (Table 1). The bismuth atom on the hollow site of the carbon hexagon has 
the smallest adsorption energy, meaning that this configuration is less stable. The bridge 
and top sites have comparable adsorption energies with the former possessing the largest 
adsorption energy. This suggests that bismuth atoms are most likely to be observed at 
bridge sites. Moreover, the optimal distance h between the adatoms and graphene surface 
is distinct for the different adsorption sites. The shorter distance at the bridge sites in- 
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dicates stronger interactions between them, and is responsible for the structural stability. 


The adatom height of 2.32 A is consistent with that of STM measurements. [50] 
Table 1; Adsorption energies and heights for various atomic sites. 


site 

AE (eV) 

h(A) 

hollow 

0.6602 

2.51 

bridge 

1.3450 

2.32 

top 

1.2904 

2.34 


The distribution of bismuth atoms can be clearly explored by calculating the ground- 
state energy. The bridge site is a much more stable one within a single hexagonal ring, 
according to the above-mentioned results. The ground-state energies for the bridge sites in 
different hexagons along the periodic armchair direction are all calculated, these energies 
provide the full information needed to determine the most stable position in monolayer 
graphene, as shown in Fig. 2. These positions can be further divided into the three regions, 
red, yellow, and gray ones. The red hexagonal region is closest to the buffer layer with a 
distance of 3.21 A; it corresponds to the lowest ground-state energy among all bridge sites. 
This energy is set to zero in order to compare it with the other bridge sites. A bit away 
from the red region, the bridge sites between red and yellow sticks have higher ground-state 
energies. The energy differences is about 17 — 23 meV. The other parts associated with the 
gray hexagons, which possesses ground-state energies comparable to each other, exhibit the 


highest energy difference in the range of 48-52 meV. The energy variations clearly illustrate 





that the bismuth atoms are hardly transported from the red to the other regions at room 


temperature. The Bi atoms are most stable at the red hexagonal rings, since at these sites 
the Van der Walls interactions with the crests of the buffer layer are strongest. The energy 
barrier of ~ 50 meV, which creates a potential well for the transportation of Bi adatoms, 
plays an important role in the dramatic change of adatom distribution during the variation 
of temperature. 

The adsorption and ground-state energies are useful to comprehend the large-scale pat¬ 
tern of Bi adatoms identihed by experimental measurements. The periodic hexagon image 
of Bi-adsorbed monolayer graphene is shown in Fig. 3(a), where the red and grey regions 
represent the most stable and unstable positions, respectively. The shortest, longest, and 
average distances between two Bi adatoms are, respectively, 14.2 A, 18.1 A, and 15.9 A. The 
latter two are indicated by brown and black arrows between each red regions. This means 
that the most stable interatomic distance lies in the range of 14.2-18.1 A. The simulated 
pattern deserves closer examination with the STM measurements, [501 [CT| which reveals 
the distribution of the Bi adatoms. The large-scale hexagonal array of Bi atoms is clearly 
identihed in Fig. 3(b). Via the statistical analysis, it is determined that the interatomic 
distance is 15 A and 16 A for most of Bi adatoms, but 14 A and 17 A for some others (inset 
in Fig. 3(b)). All the simulation results agree with these STM measurements. This further 
illustrates that the periodic corrugation of the buffer layer indeed inhuences the monolayer 
graphene, and thus the arrangement of the Bi adatoms. The stable uniform distribution in 
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room temperature should be very useful in the future applications in energy engineering. 


Besides the most stable structure, there are metastable nano-structures induced by the 
annealing treatment. [SI] As the temperature of Bi-adsorbed graphene is raised from 300 
K to 500 K and then reduces to the original one, During the treatment, many Bi-adatom 
nanoclusters are observed in Fig. 4(a). Most of these are distributed in triangular and 
rectangular arrangements. These special conhgurations can be understood by a look at the 
ground-state energies and Bi-Bi interaction energies of optimized structures, as shown in 
Table 2. The Bi adatoms are set at the bridge sites of a hexagonal ring to get stronger 
interactions, as indicated in Fig. 4(b) for various adatom numbers (from 1 to 6). The 
higher the number of adatoms is, the lower Etotai- The stronger attractive Bi-Bi interactions, 
compared to the Bi-C bonds, are responsible for this result. The latter, Esi-c = —1-01 eV, 
is obtained by subtracting the ETotai of graphene from that of single-Bi-adsorbed graphene. 
The reduced energy due to the Bi-Bi interactions can be calculated by 

n 

^EBi-Bi = {Erotai — Ecra — Esuf “ nEsi — EBi-c)/n. ( 2 ) 

i=l 

Among the various Bi-adatom numbers, AEBi-Bi^s in the 3-, 4-, 5-, and 6-adatom nan¬ 
oclusters are much lower than that of the 2-adatom one. This indicates that the former 
four are metastable structures. The experimental measurements show that most of the 
nanoclusters are composed of 3 and 4 adatoms (green and yellow arrows in Fig. 4(a)). If 
the bismuth coverage is insufficient, the 5- and 6-adatom nanoclusters are absent. A critical 
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Table 2: Total energies and Bi-Bi interaction energies for various Bi-nanoclusters. 


# Bi atoms 

Erotai (eV) 

^Esi-Bi (eV) 

1 

-2317.86 

X 

2 

-2322.14 

-1.13 

3 

-2326.77 

-1.62 

4 

-2330.60 

-1.67 

5 

-2334.19 

-1.65 

6 

-2337.86 

-1.65 


factor is that the bismuth atoms from the large-scale hexagonal array need to overcome 
the energy barrier (~ 50 meV in Fig. 2) in order to form these patterns. The temperature 
increase causes the transport of Bi adatoms between two neighboring unit cells (red hexag¬ 
onal in Fig. 3(a)), where they have the same probability of moving toward or away from the 
hexagonal unit cell. The hexagonal symmetry leads to the 4-adatom-dominated nanoclus¬ 
ters. However, the relatively few vacancies in the large-scale hexagonal array (arrows in 
Fig. 3(b)) can create a the non-uniform transport environment and thus the 3-adatom nan¬ 
oclusters. It should also be noted that the nearest distance between two bismuth clusters 
is about 16 A, revealing that the energetic favorable adsorption sites of them correspond 
to the red hexagonal rings in Fig. 2. This further illustrates that the buffer layer plays an 
important role at various temperatures. 

The density of states, shown in Figs. 5(a) and 5(b), can directly reflect the primary 
electronic properties. For a hexagonal array of Bi-adsorbed graphene, DOS is finite at 
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E = 0; it possesses a dip at low energy, and a peak structure at ~ —0.6 eV (Fig. 5(a)). 


The first feature means that there exists a certain amount of free carriers. The second one 
at —0.2 eV is ascribed to the Dirac point. The third one comes from the contribution of 
bismuth atoms. Such electronic states are critical characteristics to identify the existence 
of bismuth atoms. The STS measurements, in which the tunneling differential conductance 
map of the dl/dV-V curve is proportional to the DOS, can provide an accurate and efficient 
way to examine theoretical calculation. Previous experiments [5011^ show features similar 
to the DFT calculation, as shown in Fig. 5(b). The hnite DOS at i? = 0 and the small dip 
at low energy show good agreement with the theoretical results. The peak at ~ —0.7 V 
originating from the bismuth atoms is similar to the peak structure in our calculation. The 
above-mentioned comparison of electronic properties might promote potential application 
in electronic devices. 

4. Conclusion 

The geometric and electronic properties of bismuth adsorbed on monolayer graphene are 
studied by ab initio density functional theory calculations. All the critical interactions are 
calculated, namely the six-layered substrate, periodically corrugated buffer layer, slightly 
deformed monolayer graphene, and adatom arrangements. The ground-state energies, bis¬ 
muth adsorption energies, and Bi-Bi interaction energies for various Bi adatom conhgu- 
rations are investigated in detail. Remarkably, the corrugated buffer layer has periodic 
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non-uniform Van der Walls interactions with the monolayer graphene. Adatom distribu¬ 


tions are shown to be enriched by the substrate and the rippled buffer layer. The optimized 
structure, with an inter-adatom distance of 14 ~ 17 A, presents a large-scale hexagonal 
array. The energy barrier for the formation of the adatom nanoclusters is about 50 meV. 
With increasing temperature, 3- and 4-adatom nanoclusters are, respectively, revealed in 
triangular and rectangular patterns, being associated with the large-scale hexagonal sym¬ 
metry, Bi coverage, and sparsity of vacancies. The stable and metastable structures are 
consistent with STM measurements of the annealing processes. Furthermore, DOS exhibits 
a hnite value at the Fermi level, a dip structure at ~ —0.2 eV, and peak at ~ —0.6 
eV, as identihed by the STS measurement. These structures arise from the free conduction 
electrons, the shifted Dirac-cone, and the bismuth-dominated electronic states, respectively. 
These features are consistent with those observed by the STS measurements. The stable 
conhguration in room temperature, tunable adatom arrangements, and feature-rich elec¬ 
tronic properties may be potentially important for the application of bismuth-adsorbed 
systems in nanoelectronic devices and energy materials. 
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FIGURE CAPTIONS 


Figure 1 Geometric structure of silicon carbide substrate, buffer layer, monolayer rgaphene, 
and bismuth adatoms. C, Si and Bi are, respectively, indicated by the gray, yellow and 
purple solid circles. 

Figure 2 Ground-state energies of bismuth adsorption on different sites above monolayer 
graphene. 

Figure 3 Geometric structures of hexagonal Bi array by the (a) DFT calculations and 
(b) STM measurements. 

Figure 4 Geometric structures of Bi nanoclusters by the (a) STM measurements and (b) 
DFT calculations. 3- and 4-adatoms are, respectively, indicated by the green and yellow 

arrows. 

Figure 5 Density of states for hexagonal Bi array by the (a) DFT calculations and (b) 
STS measurements. 
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Figure 1: Geometric structure of silicon carbide substrate, buffer layer, monolayer rgaphene, 
and bismuth adatoms. C, Si and Bi are, respectively, indicated by the gray, yellow and 
purple solid circles. 
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Figure 2: Ground-state energies of bismuth adsorption on different sites above monolayer 
graphene. 
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Figure 3: Geometric structures of hexagonal Bi array by the (a) DFT calculations and (b) 
STM measurements. 
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Figure 4: Geometric structures of Bi nanoclusters by the (a) STM measurements and (b) 

DFT calculations. 3- and 4-adatoms are, respectively, indicated by the green and yellow 
arrows. 
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Figure 5: Density of states for hexagonal Bi array by the (a) DFT calculations and (b) 
STS measurements. 
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